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Abstract— Based on currently available data on the three-dimensional field of space ve- 
locities of young (<50 Myr) open star clusters and the radial velocities of HI clouds and 
star-forming (HII) regions, we have found the Galactic rotation curve in the range of Galac- 
tocentric distances 3 kpc < i? < 12 kpc using the first six terms of the Taylor expansion of 
the angular velocity of Galactic rotation in Bottlinger's equations. The Taylor terms found 
at the Galactocentric distance of the Sun _Ro = 7.5 kpc are: ojq = —27.7 ±0.6 km s~^ kpc~^, 
col = 4-13 ± 0.07 km s'^ kpc'^, tug = -0.912 ± 0.065 km s'^ kpc'^, to^ = 0.277 ± 0.036 km 
s-^ kpc"^, = -0.265±0.034 km s"^ kpc^^ = 0.104±0.020 km s"^ kpc^^ In this case, 
the Oort constants are A = 15.5 ±0.3 km s^^ kpc~^ and B = —12.2 ±0.7 km s~^ kpc^^. We 
have established that the centroid of the sample moves relative to the local standard of rest 
along the Galactic Y axis with a velocity of —6.2 ±0.8 km s~^. A Fourier spectral analysis of 
the velocity residuals from the derived rotation curve attributable to density waves reveals 
three dominant peaks with wavelengths of 2.5, 1.4, and 0.9 kpc and amplitudes of 4.7, 2.6, 
and 3.6 km s~^, respectively. These have allowed us to estimate the distances between the 
density wave peaks, 1.9, 2.4, and 3.2 kpc as R increases, in agreement with the description of 
the density weave as a logarithmic spiral. The amplitude of the density wave perturbations 
is largest in the inner part of the Galaxy, ^ 9 km s~^, and decreases to ~ 1 km s~^ in 
its outer part. A spectral analysis of the radial velocities of young open star clusters has 
confirmed the presence of periodic perturbations with an amplitude of 5.9 ± 1.1 km s~^ and 
a wavelength A = 1.7 ± 0.5 kpc. It shows that the phase of the Sun in the density wave 
is close to —tt/2 and the Sun is located in the interarm space near the outer edge of the 
Carina-Sagittarius arm. 

INTRODUCTION 

Constructing the Galactic rotation curve is of great importance in solving a number of 
problems, such as estimating the mass of the Galaxy, determining the distribution of matter, 
estimating the hidden mass, studying the dynamics of the Galaxy and its subsystems, etc. 
The shape of the rotation curve in the outer Galaxy is of particular interest. Xue et al. 
(2008) showed that the circular rotation velocity is close to 200 km s~^ far from the Galactic 
center, being 175 km s~^ at a distance of 60 kpc. 
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The following data on various objects in the Galactic disk are used to determine the 
Galactic rotation parameters: HI and HII radial velocities (Burton 1971; Clemens 1985; 
Fich et al. 1989; Brand and Blitz 1993; Nikiforov 1999; Russeil 2003; Avedisova 2005), data 
on open star clusters (OSCs) and associations, Cepheids and other young stars (Mishurov 
and Zenina 1999; Rastorguev et al. 1999; Dambis et al. 2001; Zabolotskikh et al. 2002; Dias 
and Lepine 2005). 

The most comprehensive information about the pattern of Galactic rotation in the inner 
region can be obtained using the radial velocities of hydrogen clouds extracted from 21-cm 
or CO radio observations. 

Once the theory of spiral density waves has been developed (Lin and Shu 1964; Lin 
et al. 1969), it turned out that the small deviations of the circular velocities of hydrogen 
clouds from the smooth rotation curve constructed from radial velocities could be successfully 
interpreted in terms of the density wave theory (Burton 1971; Burton and Gordon 1978; 
Clemens 1985). However, the hydrogen radial velocity data arc insufficient to trace the 
kinematic effect of the spiral pattern in the outer Galaxy. (Indeed, since the tangential point 
method is inapplicable, the hydrogen data have large errors.) Therefore, other data, for 
example, those on OSCs, should be invoked additionally. 

When only the radial velocities are analyzed, the following two quantities serve as model 
ones for constructing the Galactic rotation curve: the Galactocentric distance of the Sun 
Rq and the linear circular velocity of the Sun Vq (Brand and Blitz 1993; Russeil 2003; 
Avedisova 2005). Including OSCs for which reliable estimates of their proper motions are 
available allows us either to reduce the number of model quantities to one (Rq), because 
in this case the angular velocity of solar rotation cuq can be determined from observations 
(Rastorguev et al. 1999; Zabolotskikh et al. 2002), or to determine both Rq and Vq directly 
from observational data (Bobylev et al. 2007). 

The goal of this paper is to solve two problems: (1) constructing a smooth rotation 
curve from young objects in the Galactic disk and (2) analyzing the high-frequency periodic 
circular velocity residuals of the objects used from the derived rotation curve attributable 
to the effect of spiral density waves. 

To solve the first problem, the R range is usually broken down into several segments: 
the central (i? < 0.8 — 1.5Rq), inner (pa 1.5i?o < R < Rq), and outer {R > Rq) regions 
of the Galaxy, as was done by Burton (1971) and Clemens (1985). This approach is used 
mainly because there is a segment of solid-body rotation in the inner Galaxy that passes 
into a fairly flat curve further out. Therefore, the entire R range is difficult to describe by a 
smooth curve. In addition, there are significant radial motions of gas at i? < 3 kpc (Bania 
1977; Liszt and Burton 1980). Several authors suggested their smoothed rotation curves for 
each of these regions (Burton and Gordon 1978; Clemens 1985). However, it is not always 
convenient to use a rotation curve defined by different expressions in different R segments. 
It is much more convenient to have one analytical expression for the rotation curve in the 
entire range. 

Here, we restrict our analysis to the range of distances R > 2.5 kpc. Zabolotskikh 
et al. (2002) found a smoothed Galactic rotation curve in a wide range of distances R 
without dividing it into segments, but this rotation curve was given in tabular form, which 
complicates its use. Our objective is to determine the Galactic rotation curve in a wide 
range of distances, R > 2.5 kpc, in analytical form. 

We use a Fourier spectral analysis to analyze the velocity residuals of the objects from 
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the derived rotation curve produced by a spiral density wave. 

Apart from the radial velocities of hydrogen clouds, we invoke data on star-forming 
regions (Russeil 2003) and new data on young OSCs (Kharchenko et al. 2007). To determine 
the projections of the OSC circular velocities, we use their total space velocities. 

DATA 

The Working Data Set 

We took the radial velocities of HI clouds obtained by the tangential point method from 
Fich et al. (1989), where the 21-cm hydrogen radio observations described in Burton and 
Gordon (1978) are presented in tabular form. We also used the CO radio observations by 
the tangential point method from Clemens (1985). 

We took the data on star-forming (H II) regions mainly from the catalog by Russeil 
(2003), which contains the radial velocities and photometric distances for 204 regions. In 
this paper, we use only 89 HII regions selected according to the following criteria: 

(i) the relative random errors in the distances to the objects do not exceed 20%; 

(ii) the regions are located outside the zone with the Galactic center-anticenter direction 
\l\ < 12°, where the chance of estimating the radial velocities is reduced sharply (Russeil 
2003); 

(iii) the regions have Galactocentric distances of more than 5 kpc (closer regions are 
rejected because of the large velocity dispersion). 

For the HII region Sharpless 294, we use a new estimate of its photometric distance, 
r = 4.8 ± 0.2 kpc, from Samal et al. (2007). For the star- forming region W30H, we 
use a new estimate of its distance, r = 1.95 ± 0.04 kpc, from Xu et al. (2006) obtained 
from the trigonometric parallax (0.512 ± 0.010 mas) measured using the VLBA. For the 
HII region Sharpless 269, we use the results of high-precision radio interferometry (Honma 
et al. 2007): the trigonometric parallax 0.189 ± 0.008 mas (r = 5.28 ± 0.23 kpc), the 
radial velocity V^sr — 19.6 km s~^, and the proper motion components {Ha cos 5, /is) j2ooo — 
(—0.422 ± 0.010, —0.121 ± 0.042) mas determined relative to an extragalactic source. Thus, 
these data allow the space velocities of the region Sharpless 269 to be determined. 

For OSCs, we use the coordinates, proper motions, radial velocities, and age estimates 
from the compilation by Bobylev et al. (2007), which is based on the catalog by Piskunov et 
al. (2006). We took new radial velocities for a number of OSCs from the CRVOCA catalog 
(Kharchenko et al. 2007). In this paper, we consider only the OSCs whose ages do not 
exceed 50 Myr and for which the space velocities can be determined. Wc took the data for 
two young clusters, Cr 272 and Tr 28, which are missing from the CRVOCA catalog, from 
the catalog by Dias et al. (2002). 

Heliocentric Radial Velocities 

We will work with the radial velocity estimates for the objects given in the heliocentric 
frame of reference. Different authors give their radial velocity estimates relative to the local 
standard of rest using differing parameters — {U,V, W) lsr- 

We reduce the hydrogen radial velocities from Burton and Gordon (1978) and 
Clemens (1985) to the heliocentric frame of reference using the velocities {U,V,W)lsr — 
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(10.3, 15.3, 7.7) km s ^ . We calculated these values based on the following parameters of the 
standard solar motion: (q;,(5)i9oo = (270°, +30°) and V — km s~^ (Burton and Gordon 
1978). 

We reduce the radial velocities of the star-forming regions to the heliocentric frame of 
references using the velocities (f/, V, W)lsr = (10.4, 14.8, 7.3) km s~^ from Russeil (2003). 

We reduce the radial velocities of the Hll region Sharpless 269 to the heliocentric frame 
of reference using the velocities ([/, V, W)lsr = (10.0, 15.4, 7.8) km from Honma et al. 
(2007). 

All proper motions of the objects are given for the epoch J2000.0, i.e., they are homoge- 
neous. Therefore, no problem with their reduction to a single frame of reference arises. 

METHODS 

In this paper, we use a rectangular Galactic coordinate system with the axes directed away 
from the observer toward the Galactic center (I — 0°,b — 0°, the X axis), along the Galactic 
rotation {I = 90°, b = 0°, the Y axis), and toward the North Galactic Pole {b = 90°, the Z 
axis). 

Determining the Rotation- Curve Parameters 

The method used here is based on Bottlinger's well-known formulas (Ogorodnikov 1958) 
with the angular velocity of Galactic rotation expanded in a series to terms of the sixth 
order of smallness in t/Rq: 

Vr — —Uq cos b cos I — (1) 
—Vq cos 6 sin Z — Wq sin b— 
-Ro sin I cos b[{R - Ro)ujI/1\+ 
+... + {R-RoYuf,/5\], 
Vi — UQsinl — Vq cos / — (2) 
-{Ro cos / -r cos h)[{R - Ro)u;^/l\+ 

+... + (i?-i?o)'c^o/5!] + 

+ruJo cos b, 

Vb = Uq cos / sin 6 + (3) 
+f sin / sin 6 — Wq cos b+ 
+Ro smlsmb[{R - Ro)ujI/1\ + 
+... + {R-Ro)'uj'j5\], 

where Vr is the radial velocity; Vi = 4.7 Ar/ii cos b and H = i.Tirfj,}, are the proper motion 
velocity components in the / and b directions, respectively (the coefficient 4.74 is the quotient 
of the number of kilometers in the astronomical unit by the number of seconds in the tropical 
year); r is the heliocentric distance of the object; the proper motion components ni cos 6 and 
Hb are in mas yr"-*^, the radial velocity Vr is in km s"^; U(7),Vq,Wq are the solar velocity 
components relative to the centroid under consideration; Rq is the Galactocentric distance 
of the Sun, which we take to be Ro — 7.5 kpc (Bobylev et al. 2007); R is the Galactocentric 
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distance of the object; R, Rq, and r are in kpc. The quantity uoq is the angular velocity 
of rotation at the distance Rq, the parameters cul, ...,uJq are the derivatives of the angular 
velocity from the first to the fifth order, respectively. The distance R can be calculated using 
the expression 

i?2 ^ ^2 ^^^2 ^ _ 2^^j. cos 6 cos / + i?o . (4) 

Note also that Eqs. (l)-(3) are written in such a way that the direction of rotation from the 
X axis to the Y axis is positive. The system of conditional equations (l)-(3) contains nine 
unknowns: Uq,Vq,Wq, ujq^uJq^uJq, uJq,uJq,ujI, which can be determined by the least-squares 
method. Equations (l)-(3) can be easily generalized to determine any necessary number of 
terms of the series in the expansion of the angular velocity of rotation uJq. The system of 
equations (l)-(3) is solved with weights of the form 

P{r,i,b) = Sq/ \[W^rair^^, 

where Pr is the weight of the equation for the radial velocity, Pi and Pfe are the corresponding 
weights for the components Vi and H, and 5*0 denotes the dispersion averaged over all 
observations, which has the meaning of the "cosmic" dispersion that we take to be 8 km s~^. 
The errors crv; and (Tv^ in the velocities Vi and H can be calculated using the formulas 



'^{Vi,v,) = 4.74r 

where (T7r/7r is taken to be 0.2, a typical error in the photometric distance for OSCs. 
Determining the Projections Vrot and Vrad 

When only the radial velocity is available for an object, we calculate the projection of the 
circular rotation velocity Vrot from the well-known formula (Burton 1971) 

Vrot — I R^O I + RVr/ {Rq sin / COS 6) . (5) 

For OSCs, we have the rectangular coordinates X, Y, Z and calculate the components of 
the observed space velocities U and V from the projections Vr, VI, and Vh (Kulikovskii 
1985). Using these components, we find two projections, Vrad directed radially away from 
the Galactic center and Kot orthogonal to it, from the relations 

Vrot = Usme+ {Vo + V) cos 6, (6) 
Vrad = -U COS e+{Vo + V) siu 9, (7) 

where Vq = |i?o^o| ^-nd the position angle 6 is defined as tan 6 = Y/ (Rq—X). Wc assume that 
the velocities U and V are free from the solar velocity relative to the centroid Vq{uq, Vq, Wq) 
derived above from Eqs. (l)-(3). The errors in the projections Vrad and Vrot can be estimated 
from the relations 

(^"^{Vrot) = <^u si'^^ + CTy cos^ 9, 
(^^(Vrad) = o-f; cos^ + ay siu^ 9, 

where the errors in the space velocities U and V are denoted by au and ay- 



Analysis of the Residual Velocities 

At the next step, we consider the rotation velocities AV that are the residual ones with 
respect to the Galactic rotation curve found at the first step. 

A spectral analysis consists in applying the direct Fourier transform to the sequence of 
residual velocities in the following way: 

^(A.) = ^E^^'exp(-j-i?,j, (8) 



where AV(Xk) is the value of harmonic k of the Fourier transform, M is the number of 
measurements of the residual velocities AV^ with coordinates i = 1,2, M, and is 
the wavelength in kpc; the latter is equal to D/k, where D is the period of the original 
sequence in kpc. 

Constraints 

The constraints listed below reduce significantly the error of a unit weight when the system 
of equations (l)-(3) is solved by the least-squares method. 

For HII regions: (i) the error in the distance is no more than 20%; (ii) R > 5 kpc; (iii) 
location outside the zone with the Galactic center. anticenter direction |/| < 12°. 

For OSCs: (i) the age is younger than 50 Myr; (ii) the error in the modulus of the space 
velocity is ^ct|^ + + crf^ < 15 km s^^ (iii) Vpec = VTPTWTW^ < 100 km s'^ In 
particular, we also use condition (i), because the negative K-effect, i.e., the effect of radial 
motions (Bobylev et al. 2007), for young OSCs is less pronounced than that for older ones. 

For neutral hydrogen (HI): (i) R > 2.4 kpc; (ii) several radial velocities from Burton and 
Gordon (1978) are not used; these data are marked by the circles in Figs. 3a and 4. At 
the same time, when analyzing the residual velocities, we rejected only two points close to 
the Sun. Condition (i) is used, because a transition to the zone of solid-body rotation is 
observed on the rotation curve at i? ~ 2.4 kpc (Burton 1971; Clemens 1985). 

The system of conditional equations (l)-(3) is solved twice with an analysis of the residu- 
als using the 3a criterion (we clearly see from Fig. 4 that the points rejected by this criterion 
have large residuals). 

RESULTS 

At a fixed Galactocentric distance of the Sun {Rq = 7.5 kpc), by solving the system 
of 777 equations, we found the parameters of the solar velocity relative to the centroid, 
{uq, Vq, Wq) = (8.41, 11.41, 8.01) ± (0.49, 0.50, 0.64) km s~^ and the parameters of the an- 
gular velocity of Galactic rotation, 

ujo = -27.67 ± 0.61 km s-^ kpc'^, (9) 
00^ = +4.132 ± 0.072 km s'^ kpc'^, 
cul = -0.912 ± 0.065 km s'^ kpc-^ 
uj^ = +0.277 ± 0.036 km s"^ kpc-^ 
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Lol = -0.265 ± 0.034 km s'^ kpc-^ 
= +0.104 ± 0.020 km s'^ kpc"^. 

As a result, taking into account all constraints, we used the data for 140 young OSCs 
(including S269) with a mean age of 18 Myr, 89 radial velocities of the HII regions (Russeil 
2003), 135 radial velocities for the tangential points (CO) from Clemens (1985), and 133 
radial velocities for the tangential points (H I) from Burton and Gordon (1978). The error 
of a unit weight is ctq = 8.0 km s~^. 

As a result, the Oort constants A = 0.5/?o^d) B = uq + 0.5i?o<^0' ^ ~ 15.50 ± 0.27 
km s~^ kpc^^ and B = —12.17 ± 0.67 km s^^ kpc~^. 

Based on parameters (9), we estimated the circular rotation velocity of the solar neigh- 
borhood to be Vo = l-Ro'^ol ~ 208 ± 5 km and the period of its revolution around the 
Galactic center to be T = 27r/(7 a;o) = 221 Myr, where the coefficient 7 = 1.023 x 10"^ [km 
s-^ kpc-^]/[rad yr'^] (Murray 1986). 

Figure 1 shows the coordinates of the objects being analyzed in projection onto the 
Galactic XY plane. 

Figure 2 shows the Galactic rotation curve constructed from parameters (9); the OSC 
circular velocity projections were calculated only from radial velocities using Eq. (5). As 
can be seen from Eq. (5). the division by sin/ is made in the second term. Hence, the 
projections of the objects with / close to zero have large deviations and large errors. In Fig. 
2, the open squares mark 17 OSCs with |/| < 12°. Therefore, Eq. (5) is not used to analyze 
the OSC velocities. 

Figures 3a and 3b show the OSC circular velocity projections calculated using Eq. (6). 
The thick solid line in these figures indicates the derived rotation curve. Figures 3a and 
3b differ by the age of the OSCs used, whose velocities are denoted by the circles. Figure 
3a pertains to young OSCs (younger than 50 Myr), while Fig. 3b pertains to old OSCs 
(older than 50 Myr). We see from the figure that the velocities of both young and old OSCs 
fall nicely on the derived Galactic rotation curve; in contrast to the old OSCs, the young 
OSCs exhibit a distinct high-frequency component related to the effect of density waves. 
The dotted line in Fig. 3a indicates the Galactic rotation curve from Zabolotskikh et al. 
(2002), which was found using seven terms in the expansion of the angular velocity. From 
the figure, we see good agreement between the two rotation curves for R > 3.5 kpc. 

We solved Eqs. (l)-(3) using various numbers of terms in the expansion of cuq. A direct 
construction of the rotation curves shows that adding each succeeding term extends signifi- 
cantly the boundaries of the confidence interval. However, it can be said with confidence that 
three expansion terms are clearly not enough to construct the Galactic rotation curve in the 
range of distances 3 kpc< i? < 12 kpc (Bobylev et al. 2007), while seven terms are too many. 
For example, when the seventh term is added to Eqs. (l)-(3), terms that do not differ sig- 
nificantly from zero appear: lJq = -|-0.107± 0.086 km s~^ kpc"'' and oOq = —0.076 ±0.035 km 
s^^ kpc^'', with the shape of the curve remaining virtually unchanged. In this sense, we 
found an optimal solution. 

We can see from our comparison of the OSC circular velocities in Figs. 2 and 3a that the 
projections calculated using Eq. (6) have a smaller residual velocity dispersion with respect 
to the derived rotation curve. 

The dotted line in Fig. 3b denotes the rotation curve as derived by Dias and Lepine 
(2005) for Rq = 7.5 kpc; we added a constant term of 18 km s~^ to Vq = 190 km s~^ to 
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reconcile it with our data a.t R = Rq{Vo = 208 km s~^). 

Figure 4 shows the angular velocity of Galactic rotation determined using parameters 
(9). The circles mark the HI data points rejected by the 3a criterion when seeking for the 
rotation curve. Note, however, that these points were not rejected at the second step during 
our spectral analysis of the high-frequency deviations. Based on Fig. 4, we choose the range 
of distances 3 kpc < i? < 12 kpc as the most significant one. 

The R distribution of the residual velocities obtained by subtracting the derived Galactic 
rotation curve from the original distribution of circular velocities is shown in Fig. 5. This 
distribution corresponds to the case where no OSCs are used, i.e., here we essentially re- 
peated the analysis by Clemens (1985). The result obtained is discussed in the next section. 
Figure 5b corresponds to the entire set of objects under consideration. The residual velocity 
dispersion is 8.8 km s~^. The maximum residual velocity is 45 km s"^. As we see from the 
figure, a periodic structure of the velocity distribution attributable to the effect of Galactic 
spiral density waves is traceable. Therefore, it is quite natural to use a spectral analysis to 
investigate the residual velocities in an effort to obtain quantitative characteristics. 

We can see from our comparison of Figs. 5a and 5b that the fitting residual velocity curve 
in Fig. 5b has more distinct peaks; the similarity of the two periods in the inner Galaxy is 
better seen in this curve than in the curve of Fig. 5a. 

Figure 6 shows the power spectra of the residual velocities as a function of the natural 
logarithm of the wavelength-to-i?o ratio (lnA/i?o)- Figure 6a pertains to the residual veloc- 
ities AVi obtained by subtracting Vrot — Vo (where we found the velocity Vq = 208 km s~^ 
above), representing the deviation from a flat rotation curve, while Fig. 6b pertains to the 
residual velocities AVu obtained by subtracting the derived Galactic rotation curve from 
the original distribution of Vrot (Fig- 5b). 

The power spectra are fairly wide, but we give them only up to In A/i?o = —3.5. The am- 
plitude of the spectrum decreases rapidly with increasing frequency (decreasing wavelength). 

We can see from our comparison of Figs. 6a and 6b that taking into account the derived 
Galactic rotation curve with parameters (9) allows us to get rid of the low-frequency spectral 
component almost completely. 

Three dominant peaks with wavelengths of 2.5, 1.4, and 0.9 kpc {\n\/Ro ~ 
— 1.2, —1.8, —2.2) and amplitudes of 4.7, 2.6, and 3.6 km s~^, respectively, can be distin- 
guished in the power spectrum of the residual velocities shown in Fig. 6b. The thick line in 
Fig. 5b indicates the fitting curve of the residuals constructed using the frequencies of the 
spectrum in Fig. 6b for IhX/Rq > —2.5. In contrast to the approach by Clemens (1985), 
who estimated the perturbation wavelength and amplitude only from one dominant peak in 
the residual power spectrum, we construct the fitting curve in a wider frequency range. As 
a result, we obtain a complex periodic curve that represents satisfactorily the residuals. 

Four main peaks with the following characteristics can be distinguished in Fig. 5b: 

(1) Ri = 3.9 kpc, AVi = 8.2 km s'^; 

(2) i?2 = 5.8 kpc, AV2 = 9.2 km s"^; 

(3) i?3 = 8.2 kpc, AV3 = 1.0 km s'^; 

(4) i?4 = 11.4 kpc, AV4 = 1.8 km s-^ 

The distance differences between the peaks are Ai?2i = 1.9 kpc, Ai?32 = 2.4 kpc and 
AR43 — 3.2 kpc. As can be seen from Fig. 5b, the fourth peak is revealed with the lowest 
confidence. In general. Fig. 5b gives only a qualitative picture of the effect of spiral density 
waves. This is because the line of sight passes through different sections of the spiral arms 
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and, hence, the resulting velocity distribution is "smeared" . Nevertheless, we can suggest the 
following identification with known spiral arms: peak 1 — the Scutum-Crux arm; peak 2 — 
the Carina-Sagittarius arm; peak 3 — the Perseus arm; peak 4 — the Norma-Cygnus arm. 
In general, the locations of the first three peaks agree well with the two-arm model by Yuan 
(1969). Comparison with Figs. 4 and 5 from Russeil (2003) shows that this identification 
agrees satisfactorily with both three-arm and four-arm models of the spiral pattern (see also 
Fig. 2 from Vallee (2002)). However, for the best agreement with models, for example, the 
model by Russeil (2003), we must assume that the peak with R — 2.9 kpc and AV — 1.9 km 
s~^ is significant in our Fig. 5b. In this case, it will be either the beginning of the Perseus 
arm for the three-arm model or the beginning of the Norma-Cygnus arm for the four-arm 
model. 

Figure 7 shows the residual radial velocities of young OSCs as a function of the Galacto- 
centric distance calculated using Eq. (7) . The power spectrum of the residual radial velocities 
for young OSCs obtained using Eq. (8), where the radial velocities are considered as AV, 
is shown in Fig. 8. The curve of residuals shown in Fig. 7 has a wavelength A = 1.7 kpc 
and an amplitude of 5.9 km s~^; the linear displacement along the vertical axis is —1.6 km 
s"^. We clearly see from our comparison of Figs. 6b and 8 that, in contrast to the circular 
residual velocities, the power spectrum of the residual radial velocities of young OSCs for 
wavelengths lnA/i?o > —2.5 is simpler: in fact, there is only one distinct and symmetric 
peak near lnA/i?o ~ —1.7. Therefore, the derived curve of residuals is nearly sinusoidal in 
the range R = 6.9 kpc. 

When the circular velocities are analyzed, the error in A is ±0.3 kpc and the error in the 
perturbation velocity amphtude is ±0.5 km s"^ (504 points). We estimated these values by 
the Monte Carlo method. When the radial velocities of 140 young OSCs are analyzed, these 
errors have the following values: the error in A is ±0.5 kpc and the error in the perturbation 
velocity amplitude is ±1.1 km s~^. 

DISCUSSION 

The Galactic Rotation Curve 

Having analyzed the data on hydrogen clouds in a wide R range, up to 2Ro, Brand and 
Blitz (1993) found a nearly fiat Galactic rotation curve {A = —B). However, the Oort 
constants A = 15.50 ± 0.27 km s"^ kpc"^ and B = -12.17 ± 0.67 km s"^ kpc"^ we found, 
which are in good agreement with the result of various recent studies (for an overview of 
the determinations of the Oort constants, see Bobylev 2004), show that the rotation curve 
in the solar neighborhood is nevertheless not fiat. The low-frequency peak in Fig. 6a shows 
that the amplitude of the deviations of the Galactic rotation curve from a fiat one is about 
8 km s~^ in the range 3 kpc < i? < 12 kpc. 

Note that, in general, our Galactic rotation curve in the R range under consideration 
agrees satisfactorily with the "composite" rotation curve by Clemens (1985). 

The shape of the rotation curve in the outer Galaxy depends significantly on the adopted 
value of Rq (Clemens 1985; Olhng and Merrifield 2000; Zabolotskikh et al. 2002). Our 
Galactic rotation curve is intermediate between the two rotation curves from Oiling and 
Merrifield (2000) constructed using currently available data for Rq — 7.1 and 8.5 kpc. 
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As we see from Fig. 3a, our Galactic rotation curve agrees well, within the la confidence 
interval, with one of the best (since it was constructed by carefully reconciling the distance 
scales of various samples) curves to date constructed by Zabolotskikh et al. (2002) for 
Rq = 7.5 kpc. 

As we see from Fig. 3b, our Galactic rotation curve differs significantly from the rotation 
curve that was found by Dias and Lepine (2005) using data from Clemens (1985) only in the 
distant outer part of the Galaxy (i? > 11 kpc). This difference stems from the fact that the 
observational data in the outer Galaxy are scarce, while the available data are unreliable. 
In this respect, the unique observational data for S269 (Fig. 3a) are of great importance. 

In general, it can be concluded that our Galactic rotation curve in the range 3 kpc 
< i? < 12 kpc under consideration describes well the observational data used and is in good 
agreement with the results of other authors. 

The Phase of the Sun in the Spiral Wave 

Having analyzed the velocities of HI clouds in the inner Galaxy, Clemens (1985) found the 
deviation of the motion of the local standard of rest from the circular one with a velocity 
^Vrot — +7.0 ± 1.5 km s~^ in the direction of Galactic rotation. 

In our view, this motion can also be associated with the effect of a density wave. Let us 
consider an example from Rohlfs (1977), 

(Ti = |a| cosx, 
AKot = I /el sin X 

for the region within the corotation radius in the case of a tightly wound spiral. Here, 
X = /{R-Xq) is the phase of the spiral wave, Xq is the phase of the Sun in the spiral 
wave, (Ti is the matter density, a is the density perturbation amplitude, Jr and fg are the 
amplitudes of the radial and tangential velocity perturbation components, the center of the 
spiral arm corresponds to the phase x = (for details, see the description to Fig. 21 in the 
monograph by Rohlfs (1977)). For the case of AVrot = +7.0 km s~^ under consideration, 
the phase of the solar neighborhood in the wave will be close to — 7r/2. It thus follows that 
the Sun is located at the outer edge of the arm. 

The solar velocity components Vq{uq,Vq,Wq) we found lead us to conclude that there 
is a motion of the centroid of our sample along the Galactic Y axis with a velocity A^ot — 
Vlsr ~ Vq = —6.2 ± 0.8 km s~^ relative to the local standard of rest with parameters 
(f/, V, W)lsr = (10.00, 5.25, 7.17) ± (0.36, 0.62, 0.38) km (Dehnen and Binney 1998) that 
differs significantly from zero. Note that Dehnen and Binney (1998) determined the velocity 
of the Sun relative to the local standard of rest based on the proper motions of ~ 12000 
main-sequence stars from the Hipparcos Catalogue (1997). Our previous analysis (Bobylev 
and Bajkova 2007) of the space velocities for ^ 5000 F and G dwarfs using radial velocities 
showed that {U, V, W)lsr = (8.7, 6.2, 7.2) ± (0.5, 2.2, 0.8) km s'K Thus, the reliability of the 
parameters {U, V, W)lsr is beyond doubt. 

As we see, the velocity AVrot of the centroid of our sample along the Y axis is negative. 
Hence, in this case, the Sun should be located at the inner edge of the spiral arm, having a 
phase close to 7i/2. This contradiction between the two estimates of the phase of the Sun 
was also pointed out by Zabolotskikh et al. (2002). 

in 



(10) 



It is clearly seen from Fig. 5a, where we essentially repeated the analysis of data 
from Clemens (1985), that the perturbation velocity in the solar neighborhood is positive, 
AVrot — +2 km s~^. As can be seen from Fig. 5b, including the data on OSCs in the 
analysis changes significantly the picture in the immediate neighborhood of the Sun. The 
perturbation velocity at R = Rq is AVrot = —5.5 ± 0.5 km s~^; the phase of the Sun in the 
wave should then be close to 7i/2. 

However, we again run into a significant contradiction. If the phase of the Sun is assumed 
to be close to 7r/2, then the radial velocities should decrease with increasing distance from 
the Sun in the direction of increasing R, but the reverse is true in Fig. 7. 

Note that our OSC residual velocity distributions (Figs. 5b and 7) are in good agree- 
ment with the results of the analyses of the velocity {Vr and V^) fields for OB associations 
performed by several authors: see Fig. 3 in Mel'nik et al. (2001) and Fig. 1 in Mel'nik 
(2003). 

We associate the above contradictions with the influence of the Orion arm, since a sig- 
nificant fraction of the nearby OSCs (Fig. 1) belong to this structure. In the opinion of 
several authors (Yuan 1969; Weaver 1970; Elmegreen 1980; Efremov 1997), the Orion arm 
is a branch at the outer edge of the Cygnus spiral arm (Melnik 2003) or a spur. Previously, 
Bobylev et al. (2007) found that the centroid of the OSCs belonging to the Orion arm lags 
behind the local standard of rest with a velocity of ~ —10 km s~^, while the typical lag 
velocity for the centroids of the OSC belonging to other arms is ~ —5 km s~^. 

All of the data obtained are difficult to reconcile within the framework of a simple model. 
As was shown by Sitnik and Mel'nik (1999) and Mel'nik (2003, 2006), the real picture of the 
solar neighborhood may be much more complex. 

Nevertheless, the phase of the Sun in the spiral wave can be determined from radial 
velocities. The amplitude of the OSC radial velocity perturbations is higher than that of the 
circular velocity perturbations. The radial velocities are almost free from the discrepancy 
in the velocities of the centroid and the local standard of rest (see the vertical displacement 
in Fig. 7). The amplitude of the radial velocity perturbations is essentially determined by 
the spiral density wave. Thus, for example, as was pointed out by Lin et al. (1969), the 
stars in the solar neighborhood can go away from their birthplace to distances of ~ 1.2 kpc 
over a period of ~ 10 Myr, which is attributable mainly to the rotation velocities, while the 
stars go away in the radial direction to distances that are a factor of 10 smaller due to the 
small pitch angle of the spiral arms. Based on the data in Fig. 7, we conclude that the 
phase of the Sun in the density wave is slightly larger than — 7r/2 (a more accurate value is 
Xq = — 0.667r). Consequently, the Sun is located in the interarm space near the outer edge 
of the Carina-Sagittarius arm. This location of the Sun agrees well with the distribution of 
stars and gas (Vallee 2002; Russeil 2003) and with predictions of the density- wave theory. 
Indeed, as we see from relations (10), the variations in velocity Vrad are in phase with the 
variations in density ai. 

The Scale Factor A 

The distance between the neighboring spiral arms along the Galactic radius vector is usually 
denoted by A; occasionally, it is called a scale factor or scale length (Mel'nik et al. 2001). 

Clemens (1985) obtained an estimate of A = 0.22i?o- For the value of Rq = 8.5 kpc used 
by him, this gives A = 1.9 kpc, which is in good agreement with our results. 
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Having analyzed the space velocities of OB associations within 3 kpc of the Sun, Mel'nik 
et al. (2001) estimated the scale length of the periodic variations in radial residual velocity 
components along the Galactic radius vector to be A = 2.0 ± 0.2 kpc. 

At the same time, we found that the distance between the spiral arms does not remain 
constant in the R range under consideration and is Ai?j = 1.9,2.4,3.2 kpc. This is in 
agreement with the description of the density wave as a logarithmic spiral, as distinct from 
an Archimedean spiral characterized by a constant distance between the neighboring turns 
of the spiral. Details of the discussion on the specific shape of the spiral in the Galaxy can 
be found in the monograph by Efremov (1989). The model using an Archimedean spiral was 
suggested by Cowic and Rybicki (1982). 

Having analyzed the radial velocities of young OSCs, we obtained an independent esti- 
mate of A = 1.7 ± 0.5 kpc from data in the range of distances 6 kpc< R < 9 kpc. 

The Perturbation Amplitudes 

The perturbation amplitudes of the circular, fg = 4.6 ± 0.5 km s~^ (the main peak in Fig. 
6b), and radial, fu — 5.9± 1.1 km s~^ (the main peak in Fig. 8), velocities that we found are 
in good agreement with the minimum perturbation amplitude of 5 km s~^ (at the minimum 
wavelength A = 0.22i?o) estimated by Clemens (1985). 

Based on formulas from Lin et al. (1069), Burton (1971) estimated the change in the 
amplitudes of the perturbation velocities from a density wave {fn and f^) as a function of R 
in the range 0.2Ro < R < 1.3it!o- It turned out that for both velocities, there is a wide peak 
in the region 0.5it!o < R < 0.9it!o, where the velocities reach ^ 8 km s~^ {fji is larger than 
/g by ~ 1 km s~^ everywhere), the velocities are about 4 km s~^ near 0.2i?o and decrease 
to 2.5 km s~^ near 1.3Rq. Our results (Fig. 5b) agree well with the estimates by Burton 
(1971). 

However, it is hard to say that our results are in good agreement with those of other 
authors obtained for the solar neighborhood. Thus, for example, having analyzed OB asso- 
ciations, Mel'nik et al. (2001) found the velocities of the perturbations from a spiral density 
wave to he fji = 7 ± 1 km s~^ and = 2 ± 1 km s^^. Popova and Loktin (2005) obtained 
completely different velocities of the perturbations from a spiral density wave, /ij = — 4 ± 5 
km s~^ and = 13 ± 3 km s~^, using data on OSCs and OB stars. Based on Cepheids, 
Mishurov and Zenina (1999) found fn = 3.3 ± 1.6 km s"^ and fe = -7.9 ± 2.0 km s'^ The 
analysis of data on Cepheids performed by Popova (2006) showed that the amplitude of the 
perturbation velocities, J'r = — 1.8 ± 2.5 km s^^ and fe = 4.0 ± 3.4 km s~^ are small and do 
not differ significantly from zero. Zabolotskikh et al. (2002) found -7 ± 2 km s"^ and 
/() — 1 ± 2 km s~^ from data on Cepheids and OSCs. 

In general, our results lead us to conclude that the curve in Fig. 5b describes well the 
large-scale (grand design) spiral structure; there are peculiarities at R — Rq attributable to 
the kinematic properties of the Orion arm. 

CONCLUSIONS 

Based on currently available data on the three-dimensional field of space velocities of young 
(< 50 Myr) open star clusters and the radial velocities of HI clouds and star-forming (H II) 
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regions, we constructed the Galactic rotation curve in the range of Galactocentric distances 
3 kpc < i? < 12 kpc. 

At the first step, we determined a smoothed curve that fitted fairly accurately the Galactic 
rotation curve in this range of distances. The curve parameters (9) were found by using the 
first six terms of the Taylor expansion of the angular velocity of Galactic rotation at a given 
Galactocentric distance of the Sun, Rq = 7.5 kpc. We found that the centroid of the sample 
moves relative to the local standard of rest along the Galactic Y axis with a velocity of 
-6.2 ±0.8 km s'^ 

At the second step, we performed a spectral analysis of the circular velocity residuals 
of the objects under consideration from the derived rotation curve, which attributes the 
observed periodic residuals to the effect of density waves. Our spectral analysis showed that 
the peak with an amplitude of 4.6 ± 0.5 km s~^ corresponding to a wavelength A = 2.5 ± 0.3 
kpc is dominant. A similar analysis of the radial velocities for young OSCs confirmed the 
presence of periodic perturbations from a density wave with an amplitude of 5.9 ± 1.1 km 
s~^ and a wavelength A = 1.7 ± 0.5 kpc. 

Based on the described approach, Clemens (1985) obtained similar results, but only 
for the inner Galaxy (i? < Rq), where highly accurate data on hydrogen are available for 
tangential points. In contrast, for the outer Galaxy {R > Rq), the radial velocities of 
hydrogen (HI and HII) clouds usually have large errors, which makes it much more difficult 
to interpret the residuals. 

We were able to advance slightly farther than Rq, because we used the total space ve- 
locities of OSCs to determine the projections of their circular velocities. In our view, this 
approach improves the homogeneity of the sample.the tangential point method allows the 
total space velocity of hydrogen clouds to be analyzed, since the velocity vector lies entirely 
on the line of sight under the assumption of purely circular rotation of hydrogen at the 
tangential point. Therefore, supplementing the sample of objects with OSCs seems natural. 

We found that the distance between the spiral arms increases with Galactocentric dis- 
tance, being 1.9, 2.4, and 3.2 kpc; this is in agreement with the description of the density 
wave as a logarithmic spiral. The perturbation amplitude also changes — it is largest in the 
inner part of the Galaxy. ~ 9 km s^^, and decreases to 1 km s^^ in its outer part. 

Analysis of the radial velocities of a sample of young OSCs suggests that the phase of 
the Sun in the density wave is close to — 7r/2 and that the Sun is located in the interarm 
space near the outer edge of the Carina.Sagittarius arm. 
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Fig. 1. Coordinates of the objects being analyzed in projection onto the Galactic XY plane. 
The Sun lies at the coordinate origin. The pluses, dots, and filled circles denote the HII 
regions, tangential points, and OSCs, respectively. 
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Fig. 2. Galactic rotation curve. The vertical line indicates Rq = 7.5 kpc. For OSCs, the 
velocities were calculated only from the radial velocities using Eq. (5). 




Fig. 3. Galactic rotation curve (thick line), the thin hnes mark the boundaries of the la- 
confidence intervals, the vertical line indicates Rq = 7.5 kpc, the symbols in panel (a) 

mark the tangential points that were not used to determine the rotation curve parameters 
and the remaining symbols are the same as those in Fig. 1. The dotted line in panels (a) 
and (b) denote the rotation curve as derived by Zabolotskikh et al. (2002) and Dias and 
Lepine (2005), respectively. 
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Fig. 4. Angular velocity of Galactic rotation versus Galactocentric distance. The symbols are 
the same as those in Figs. 1 and 3, the thin lines mark the boundaries of the la confidence 
intervals, and the vertical line indicates Rq — 7.5 kpc. 
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Fig. 5. Residual circular velocities versus Galactocentric distance: (a) no OSCs are used; 
(b) the entire set of objects being analyzed is used. 
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Fig. 6. Power spectrum of the residual rotation velocities as a function of ln(A/i?o) 
residual velocities relative to the fiat curve; (b) residual velocities relative to curve (9) 
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Fig. 7. Residual radial velocities of young OSCs versus Galactocentric distance. 




Fig. 8. Power spectrum of the residual radial velocities for young OSCs as a function of 
ln(A/i?o). 
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